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A model of multi-component hybrid inflation, dubbed multi-brid inflation, in which various ob- 
servable quantities including the non-Gaussianity parameter /jvl can be analytically calculated was 
proposed recently. In particular, for a two-brid inflation model with an exponential potential and the 
condition that the end of inflation is an ellipse in the field space, it was found that, while keeping the 
other observational quantities within the range consistent with observations, large non-Gaussianity 
is possible for certain inflationary trajectories, provided that the ratio of the two masses is large. 
One might question whether the resulting large non-Gaussianity is speciflc to this particular form 
of the potential and the condition for the end of inflation. In this paper, we consider a model of 
multi-brid inflation with a potential given by an exponential function of terms quadratic in the 
^^ ' scalar field components. We also consider a more general class of ellipses for the end of infiation 

*vj I than those studied previously. Then, focusing on the case of two-brid inflation, we find that large 

non-Gaussianity is possible in the present model even for the equal-mass case. Then by tuning 
r^^ • the model parameters, we find that there exist models for which both the non-Gaussianity and the 

P~r' ' tensor-to-scalar ratio are large enough to be detected in the very near future. 
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I. INTRODUCTION 



The primordial non-Gaussianity has been one of the hottest topics in cosmology in recent years. The conventional, 
single- field slow- roll inflation predicts that the curvature perturbation is Gaussian to an extremely high accuracy [1]. 
In other words, if any primordial non-Gaussianity is detected, it strongly indicates that the dynamics of inflation is 
^ I not as simple as we thought it to be. 

The primordial non-Gaussianity is conveniently represented by a parameter denoted by f^^ [2]. Roughly, it is the 
-- -. ' ratio of the 3-point correlation function (or the bispectrum) to the square of the 2-point correlation function (or the 
►^ . square of the spectrum). It is expected that near- future experiments such as those of PLANCK will be able to detect 
^^ ' Jnl at a level as small as 5 [3]. 

00 . Finding even a small deviation from Gaussianity will have profound implications on the theory of the early universe. 

T-H ■ Consequently, numerous types of inflationary models that produce detectable non-Gaussianity have been proposed and 
^D [ studied [4-12]. In terms of the nature of non-Gaussianities, most of these models can be classified into two categories; 
[~^ • those with non-Gaussianities arising intrinsically from the quantum fluctuations, and those with non-Gaussianities due 
f^ ' to nontrivial classical dynamics on superhorizon scales. A typical example of the former is the DBI inflation, in which 
25 , the slow-roll condition can be fully violated [13]. In this case, the equilateral /nl (denoted by /^^' ) representing 
the amplitude of the bispectrum of the equilateral configurations, is found to play an important role [14]. On the 
other hand, in the latter case where non-Gaussianities are produced on superhorizon scales, by causality the local 
/nl (denoted by /]vl*') characterizes the level of the non-Gaussianity. It is defined in terms of the coefficient in front 
of the second order curvature perturbation [2] , 



* = *L + ./]^r'*i , (1-1) 

where $ is the curvature perturbation on the Newtonian slice and ^l is its linear, Gaussian part. 

In this paper, we focus on the latter case, that is, we consider models that may produce a large value of fj^j^^, 
for example 10-100. More specifically, we consider hybrid infiation with multiple infiaton fields, dubbed multi-brid 
infiation [15]. The infiaton fields are assumed to follow the slow-roll equations of motion, and their fluctuations are 
assumed to be Gaussian. In this case, the SN formalism is most useful for the evaluation of the curvature perturbation 
and non-Gaussianity [16-20]. 

As in the conventional hybrid inflation, the inflaton fields are coupled to a water-fall field, and inflation ends when 
the infiaton fields satisfy a certain condition that triggers the instability of the water-fall field. However, unlike the 
case of a single inflaton field in which there is essentially no degree of freedom in the condition for the end of inflation, 
there is a substantial increase in the degree of freedom at the end of inflation in multi-brid infiation and it widens the 
viable range of the parameter space considerably and leads to the possibility of generating large non-Gaussianity. 

As a model of multi-brid infiation, an analytically solvable two-brid infiation model was recently investigated in 
detail [15], where the potential was assumed to be exponential with the exponent given by a linear combination of the 
infiaton fields. In this paper, we consider a two-brid model with again an exponential potential but with the exponent 



given by a quadratic function of the inflaton fields. The potential has point symmetry about the origin of the field 
space, in contrast to the case of the linear exponent which has no symmetry. Thus by investigating the quadratic 
case, we will be able to see if the generation of large non-Gaussianity is a generic feature of multi-brid inflation or if 
it is due to the lack of symmetry that leads to large non-Gaussianity in the linear exponent case. 

In passing, we mention that the possibility of large non-Gaussianity from loop correction terms in the perturbative 
expansion has recently been discussed by CogoUo et al. [21]. Although such a case certainly needs further investigation, 
in this paper we concentrate on the case in which leading order (tree) terms dominate over loop correction terms. 

This paper is organized as follows. In §2, we describe our model and derive basic formulas to be used in the 
proceeding sections. In §3, using the formulas derived in §2, we analytically compute the spectrum of the curvature 
perturbation Vs{k), its spectral index ng, the tensor-to-scalar ratio r, and the non-Gaussianity parameter f^'ff^- In 
§4, as a couple of special cases of our model, we analyze in detail the equal-mass case as well as the case of large 
mass ratio. We find that a large f§'ff^ is possible in both cases. In particular, in the equal-mass case, by tuning 
the parameters to some extent, we find that it is possible to have both r and f§'ff^ large enough to be detected. 
We conclude the paper in §5. Some computational details are described in Appendix A. For comparison, we also 
summarize the result for the case of the linear exponent two-brid model in Appendix B. We use the Planck units 
where M','^ = SnG = 1. 



II. TWO-BRID INFLATION WITH APPROXIMATELY QUADRATIC POTENTIAL 

We consider a two-component scalar field whose action is given by 
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df,(j)Ady(j)A + V{(l)) 



(2.1) 



where the potential is given by 



V ^ Vo(x,0i,02)exp 



{m\(l)\ +ml(j)l) 



(2.2) 



with Vq being a function of a water- fall field x as well as of 0i and (f>2 , but is assumed to be constant in time during 
inflation; see Eq. (2.12) below. The same model was previously analyzed [22]. However, the condition for the end of 
inflation considered then was not general enough to allow the possibility of large non-Gaussianity. 
The Friedmann and the fleld equations are 
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(2.3) 
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where H = a/a and the dot ' denotes a derivative with respect to the cosmic proper time; ' = d/dt. The slow-roll 
equations of motion are obtained by neglecting the kinetic term in the Friedmann equation and the second time 
derivative in the field equations. Thus the slow-roll equations of motion are 
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where the number of e-folds counted backwards in time, dN = —Hdt, is used as the time variable for later convenience. 
Note that the effective mass squared M^ for each (pA is given by 
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(2.6) 



Thus, the slow-roll condition is satisfied if tti^ <C 1 and (j)\ are not too much larger than unity. Incidentally, since 
m\4>'\ ^ 1 under this assumption, the difference between the present potential and a pure quadratic potential, 
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(2.7) 




FIG. 1; Definitions of parameters a and 7 in fieid space. Tfie eliipse represents tlie surface of the end of inflation. 

is almost negligibly small (that is, they are equivalent to the leading order in the slow-roll approximation). An 
analytical solution for this separable potential model was first discussed by Starobinsky [23]. 
Introducing new field variables g^ as 
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the slow-roll equations become 



Hence we immediately obtain 
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where the number of e-folds is set to zero at the end of inflation and (J)aj is the final value of the inflaton fields 



We assume that inflation ends at 

a^ = G((/)i, (j)2) = 5i(</>i cosa -I- 02 sina)^ -I- gl{—(f>i sina -I- (f>2 cosa)^ , 
which is realized by the potential Vq given by 
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We parametrize the scalar fields at the end of inflation as 



namely, 
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(2.11) 



(2.12) 
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Figure 1 shows the definitions of the angles a and 7. The elHpse describes the surface at the end of inflation, defined 
by Eq. (2.11). The angle a describes the amount of rotation of the ellipse relative to the (/)i and (/)2 axes. The angle 
7 describes the position of the inflaton trajectory at the end of inflation. 
Since is a constant of motion, we have 
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This equation determines the parameter 7 in terms of 0i and (j)2'. 7 = 7(01,02). Hence, from Eq. (2.14), 0ij and 
02 J become functions of 0i and 02; 

(/>!,/= 01 j(01, (^2), (f>2J =(f>2ji(l)l,(t>2)- (2.16) 

With this understanding, the number of e-folds given by Eq. (2.10) becomes a function of (0i,02)- It is then 
straightforward to obtain 6N to full nonlinear order. It can be straightforwardly calculated as 

SN = iV(0i + 501, 02 + S(f>2) - A^(0i, 02) . (2.17) 

Before closing this section, let us make a small comment. As mentioned in [15], the above formula for SN neglects 
the fact that the surface at the end of inflation, determined by Eq. (2.11), is not an equipotential surface. This will 
give rise to an additional correction to the final SN. Nevertheless, it turns out that the correction is small and can 
be neglected, as discussed in [15]. 

III. CURVATURE PERTURBATION AND NON-GAUSSIANITY 

In this section, we compute the curvature perturbation of our model explicitly, and evaluate the curvature pertur- 
bation spectrum Vs, the spectral index ns, the tensor-to-scalar ratio r, and the non-Gaussianity parameter /j^lf''. 

We expand the SN formula Eq. (2.17) to the second order in (50 for Af(0i,02) given in Eq. (2.10). Note that Sj 
must be expressed in terms of Scj) with second-order accuracy. Details are deferred to Appendix A. The result is 
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where we have introduced the quantities and defined as 
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We note that tan/3 = 52/<?i- For example, in the case of Fig. 1, tan/3 is the ratio of the semiminor axis to the 
semimajor axis (hence (3 < 7r/4). 

We assume that the scalar field fiuctuations (50i and (502 are Gaussian with the dispersion, 

{S^^S^B)k^(^] Sab, (3.2) 



where t^ is the horizon-crossing time of the comoving wave number /c, where k ~ Ha. Then the curvature perturbation 
spectrum is given by 

' ■^^ ^ - (27r)3 "^^ ^ a^in4YZ~mlXW)^\ ^ H2^A, ' 

where Nk is the number of e-folds at the horizon crossing. Nk — N{th). Using the fact that 3H^ = V, the spectral 
index is found to be 

"^^ ^ e-2™?Nfcp^2_,_g-2mi7Vfcy2 p" siu^ /3 COS^ /? ' ^'^ 

The tensor-to-scalar ratio is given by 
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Now we evaluate the non-Gaussianity in our model. For convenience, we introduce the linear curvature perturbation 
TZl and the linear entropy perturbation S", 

XW^^ YZ^^ YZ^^ XW^^ 

7 — d0i + -— d02 -— 001 H 7 — 002 

T-, ^ 01 02 o ^ _02 01 (o fy\ 

^ mJYZ^mjXW ' m^ FZ - mf XVF ' ^ ^ ^ 

For the Gaussian fluctuations S(I)a given by Eq. (3.2), we see that S has the same spectrum as the curvature pertur- 
bation TZli but is orthogonal to it, 

(7^L•^>=0. (3.7) 

In terms of TZl and S, the nonlinear SN in Eq. (3.1) is re-expressed as 

6N^nL + Ifi^f 7^i +CRsnLS + CssS^ + --- , (3.8) 

5 

where the non-Gaussian parameter /j^^**' is given by 
/!f?' = -f 2 
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This is one of our main results. 

Before closing this section, let us briefly discuss the non-Gaussianity due to the other second order coefflcients Crs 
and Css- Although the derivation of their expressions is straightforward, we do not present their explicit expressions 
here, because they are as complicated as Eq. (3.9) and because they are unnecessary for the purpose of the present 
paper. We only mention that an inspection of the resulting expressions reveals that they can never become much 
larger than /]^£'^'. To be a bit more precise, their values can become large for certain ranges of the parameters, but 
when they become large, f^^'ff' also becomes large simultaneously. Hence, as discussed in [15], as long as we focus on 
the bispectrum, we can neglect their contribution. 

IV. CASE FOR LARGE NON-GAUSSIANITY 

Since Eq. (3.9) for /j^f,'*' is very complicated, it is not easy to study all possible cases in detail. However, there are 
some limiting cases in which we have a substantially simplifled expression for /j^|f' but which are yet sufRcicntly of 
interest. 

One case of interest is when the two masses are equal, mi = 7712. In this case, the potential during inflation is 0(2) 
symmetric. This symmetry is broken at the end of inflation because of condition (2.11), unless gi = 92- This model 
was discussed by Alabidi and Lyth [10] as a new mechanism of generating curvature perturbations. Another case of 
interest is when the ratio of the mass parameters are large, for example, mi ^ m2. In this large mass ratio limit, an 
inspection of the Eq. (3.9) suggests that a large value of f^§ff^ may be possible if the parameter W is very small. In 
this section, we investigate these two cases in detail. 



A. Equal mass 

First let us consider the equal mass case, m^ = m^ = rn'^- This means that there is 0(2) symmetry during inflation, 
and the symmetry is spontaneously broken at the end of inflation [10]. 

In the equal-mass case, the formulas derived in the previous section simplify considerably to 
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Note that the a-dependence has disappeared because of the symmetry. 

As is clear from Eq. (4.4) for fj^'ff, in order to obtain large non-Gaussianity, it is necessary for the factor in the 
curly brackets to become large, that is, cos 2/3 sin 27/(1 + cos 2/3 cos 27) ^ 1. This is possible cither in the limit 
(/3,7) -^ (0, 7r/2) or (/3,7) -^ {tt/2,0). Since these two limits are equivalent, let us take the hmit (/3,7) — > (0,7r/2). 
This corresponds to the situation in which the ellipse is highly elongated and the inflaton trajectory hits the ellipse 
close to one of the tips of the majoraxis. Then, setting tt/2 ~ j — S, wc obtain 

(4.5) 
ns-l ^ 2m'-- [1 + ^] , (4.6) 
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(4.8) 

To investigate in more detail the theoretical predictions of this model, let us derive expressions for r and f^^f^^ in 
terms of the observational data as much as possible. We first fix the amplitude of the spectrum Vs- The WMAP 
normalization [24] gives 

T^s^-(^) -2.5xl0-^ (4.9) 

at the present Hubble horizon scale. Also, the WMAP 5-year analysis [25, 26] gives the spectral index, 

"s = 0.96+jj;0;|. (4.10) 

Below we replace Vs and ns by these observed values. 
Noting the fact that 

iH' = V^Vo^^, (4.11) 

we obtain, from Eq. (4.9), the equation, 

a^ ^Qtt'^ xVs\r. (4.12) 

Inserting this into the square of Eq. (4.7), we obtain 
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Eliminating r from Eq. (4.6) by using Eq. (4.13) gives 
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Plugging this back into Eq. (4.13), we obtain the expression for r, 
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Also, plugging Eq. (4.14) into Eq. (4.8), we obtain 
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Assuming S^ 3> /32, and using the observed values given in Eqs. (4.9) and (4.10), the above expressions for r and 
/J^^*"' reduce to 

r ^ 2.7 X 10^(m2 + 0.02)2 ^ _j!_ , (4.17) 

f^r - 1.4 X 10«(m2 + 0.02) ^ -^^ . (4.18) 

Another useful expression may be obtained by combining the above two expressions: 

/--^(^)(^);;?^' ("») 

This tells us that for m^ > 0.02, in the very near future, both r and /)^£^' may be large enough to be detected upon 
tuning the model parameters to some extent. 

Ill Figs. 2 and 3, we show f]^^^^ and r, respectively, as functions of f3 for several different values of m^. The coupling 
constants are set to A/g** = 1. The spectral index is set to ns = 0.96, but we find the dependence of it is weak in the 
range 0.94 < ng < 0.98. In Fig. 2, each curve up to its peak is well approximated by Eq. (4.18). In both figures, if 
we vary A/g^, the curves will scale inversely proportional to A/5^. As we can see, although the values of f^^'j^^ and r 
are relatively sensitive to the values of m^ and /?, there indeed exist models with large /]^£^' and r simultaneously. 



B. Large mass ratio 

Here, we consider the case of large mass ratio. Let us tentatively assume that m^ 3> iri^. Then an inspection of 
Eq. (3.9) for f^j^]^^ suggests that a large f]^]^^ is possible if Vt^ ^ 1. Hence, let we set W ~ for simplicity and 
investigate this case in detail. We note that the only assumption we adopt is the condition W = 0; we do not assume 
a large mass ratio in the following analysis. Namely, the formulas derived below are valid for any mass ratio unless 
otherwise stated. 

The condition 1^ = implies the following relation between the model parameters: 

— tana tan7 = tana tan/3 tan7 = 1 . (4.20) 
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If tan/3 ~ 1, this means that the inflaton trajectory arrives at the ellipse along the 4>2 axis. In this case, Eqs. (3.3), 
(3.4), (3.5), and (3.9) respectively reduce to, 

^ / 9 y l + cos2/3cos27 /gy 
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FIG. 2: Non-Gaussian parameter //v£^' as a function of /3 for several different values of m^. The coupling constant parameters 
are set to \/g — 1. The spectral index is set to ns = 0.96. The curves are, from the one with the highest peak to that with 
the lowest peak, for m^ = 1/20, 1/23, 1/25, 1/27, and 1/30. 
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FIG. 3: Tensor-to-scalar ratio r as a functions of /3 for several different values of rv? . The other parameters are the same as in 
Fig. 2. The curves are, from the left to the right, for -xr? = 1/20, 1/23, 1/25, 1/27, and 1/30. 
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We note that the equal-mass Hmit discussed in the previous subsection can be obtained by setting m^ — m^ = m^ in 
the above equations, because the condition W = becomes irrelevant in the equal mass limit. 

Eq. (4.24) implies that we may have large non-Gaussianity if mf ^ m^ and/or cos 2/3 sin 27/(1 + cos 2/3 cos 27) ^ 1. 
We also note that in both cases the value of /]^£^' will be positive. This result is the same as the equal mass case and 
similar to the case of the linear exponent potential model discussed in [15]. We suspect that this positivity property 
may be generically true for all models that are capable of producing large local non-Gaussianity. 

First, let us assume that cos 2/3 sin 27/(1 -I- cos 2/3 cos 27) is of the order of unity. Recall that we have mf <g; 1 and 
1712 ^ 1 from the slow-roll condition. Then in order to obtain a large /]v£^', say fj^j^^ ^ 50, we need to have an 
extremely large mass ratio, ml/rn^ ^ 50?tij^^ ^ 1. Then, Eq. (4.22) implies that r must be extremely small, since 
we must have ng — 1 <^ 1. Therefore, large non-Gaussianity is possible only in models of very low energy inflation. 

In order to look for the possibility of both large r and large f^^'j^^, we consider the case of cos 2/3 sin 27/(1 -t- 
cos2/3cos27) ^ 1. As discussed in the previous subsection, this is realized either for (/3,7) — > (0, 7r/2) or (/3,7) -^ 
(7r/2,0). Again, since both limits are equivalent, we focus on the limit (/3, 7) -^ (0,7r/2). In this limit, setting 
7r/2 — 7 = (5 again, we have 

ns-l ^ iml-"- (<e^«-^l)^^ ^ + 1 V (4-26) 
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extremely small if ml ^ 7712. In other words, having a large mass ratio does not help in enlarging the parameter 
region in which both r and /Jv£^' are large 



Now, if we have /3 <C i5 <C 1, we can obtain large /jvL^'- However, again, Eq. (4.26) for ns — I implies r must be 



V. CONCLUSION 

We analytically investigated the curvature perturbation and its non-Gaussianity in a model of multi-field hybrid 
inflation, dubbed multi-brid inflation. The model we considered is a two-fleld hybrid inflation (two-brid inflation) 
model with the potential mimicking conventional quadratic potentials. The new ingredient of the model is the 
generalization of the condition for the end of inflation. We considered a very general coupling of the two inflaton 
fields to a water-fall field. 

Then, using the 6N formula, we derived an analytical expression for the curvature perturbation. Based on this 
expression, we obtained the curvature perturbation spectrum Vs, the spectral index ns, the tensor-to-scalar ratio r, 
and the non-Gaussian parameter /]^£'*'. We found that a large positive /j^^'^' is possible in this model. Then, at least 
for a certain limited range of the parameters, we explicitly showed that it is possible to have large non-Gaussianity 
while keeping the values of the other quantities consistent with those of the observation. In particular, we showed that 
when the two inflaton masses are equal, the parameters can be tuned so that they lead to a fairly large tensor-to-scalar 
ratio, r ^ 0.1, as well as a large non-Gaussian parameter, f]^^^^ ^ 50. These values will be at a detectable level in the 
very near future. On the other hand, interestingly, we found that having a large mass ratio in the present model does 
not help in producing both r and f]^'£'^ large enough to be detected. This is in contrast to the model studied in [15]. 

The standard lore has been that f]^^^^ is too small for models with large r or vice versa. We have shown, in this 
paper, not be the case, particularly in this model of spontaneously symmetry breaking at the end of inflation. This 
may be the most important conclusion of this work. At the moment, we have no clear physical explanation for this 
result. We hope we will be able to answer this question in the near future. 
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APPENDIX A: SN TO SECOND ORDER 

Here we evaluate SN to the second order in the perturbation. We assume the field fluctuations S(f)i and S(p2 are of 
linear order. 

First, we express the perturbation in the orbital parameter 7 in terms of S(j)i and (502- Setting ^7 = ^17 -I- S2J, 
where 61J and 62J are of linear and second orders, respectively, we take the perturbation of Eq. (2.15) to the second 
order. We obtain 

1 (501 1 (50f\ /I S(j)2 1 S(t)l 

mf (j)i 2ml ^1 / \"^2 4'2 27712 '?^2 

The linear part of the above equation determines Sij. We find 

1 501 1 (502 "^2 r i ™1 r i 

— 2~Z 2~r~ -;— 'J01 —002 

Here, for notational simplicity, we have introduced X , Y, Z and VF, which are defined by 

Y 9192 , 1 / . . X 

A = 01,/ = — (52 cosacos7 — f/i smasm7j , 

V 5152 d 1 . . 

y = -r— 01 f = -(02 cos(3!sm7 -f oi smacos7j , 

ga d-f ■■' g 

rj 9\92 , 1 . . , .X 

z = 02 f = — (92 sinQ;cos7 + gi cosasin7J , 

.9CT " g ' 

f) 1 

W = —'- — — 02 f = -(^2 sinQ!sin7 — gi cosQ!COS7) , (A3) 

ga d-i ' g 



where g = Vyi+sf- The factor 1/g in front of each of these quantities has been inserted for later convenience. 
Then, collecting the second-order terms in Eq. (Al), we find 

' ln<Pi,/ 2"3~2^'^'^2,/ (di7) 



1 m?0^ ' 771002 V"^l ^7^ "^2 ^7^ 

027 = o 



^m^ 07 ?7i2 07 



^2,,V fmi.^Y ( 2^^ + ^^ 2^"" + ^ 



2^ -r Kr 1 , .2 



2 25" 2^^ 

Too nil — 



(-517)^ 



(A4) 



where we note that 



X2 + r2^4cos2a-H4sin'«. Z^ + M/2 ^ 4 sin' « + 4 cos^ « . (A5) 

5 9 9 9 
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Now, we compute SN. Although it is straightforward to expand Eq. (2.10) to the second order in the field fluctu- 
ations, the calculation is much simpler if we take the perturbation of either of the solutions (f>i or 02 of the slow roll 
equations of motion (2.5). For example, the solution for i^i is expressed as 

0i = 0ij:e™i^ ^ iV=^(ln(/.i-ln(/iij) . (A6) 

The perturbation of the second equation gives 
1 



SN = , 

TTtJ 



»1 



d'-f 



In0ij(5i7-- i-^j - ^^^'^(t'^J {^ilf - ^^'^(I'lJ ^21 



(A7) 



Inserting Eqs. (A2) and (A4) into Eq. (A?), we obtain 

— ^ ^w^ + 2 ^T :^ 2 T^ w\^ ■ ^ ^ 

mi mi — mi, mi — ( ^2 ^ ry,^^^ ^ 

^X ^Z ^X ^Z \ ^X ~ ^Y 

Finally, we mention that we can divide SN into two contributions: one from during inflation up to a surface of 
constant potential energy, SN^, , and the contribution from the end of inflation, SN^ . In the case of the exponential 
potential model considered in [15], there was no non-Gaussianity in SN^, to the lowest order in the slow-roll parameters. 
In contrast, there exists non-Gaussianity in SN^ in the present model. Nevertheless, it can be easily shown that it is 
of the order of the slow-roll parameters, and hence is negligibly small. 

APPENDIX B: LINEAR EXPONENTIAL POTENTIAL MODEL 

In this appendix, we consider the case of a linear exponential potential, 

V = Vo expimicpi + m24>2) , (Bl) 

with the condition for the end of inflation given by 

o"^ = .9i(0i cos a + 02 sin a)^ -I- g2{~<Pi sin a + 02 cosa)^ . (B2) 

This model was discussed in [15]. However, it was assumed that a = 0. Here, for the sake of completeness, we consider 
the general condition adopted in the main text. 

As in §2, we parametrize the scalar field at the end of inflation as 

— cos 7 = 01 / cos a + 02 / sin a , — sin 7 = — 0i / sin a -I- 02 / cos a , (B3) 

9i ' ' 92 



or, conversely, 



6i f = ((72 cos a cos 7 — 51 sin a sin 7) , 02 f = (<?2 sin a cos 7 -I- <?i cos a sin 7) . (B4) 

5152 ' 5152 



Also, as before, we introduce g ~ \/ 9\ + 5^ , and X, y, Z and W as 

^ 51.92 , 1 / . . ^ 

A = 01 f = — (92 cosacos7 — oi smQ;sm7) , 

9<^ 9 ' 

V- 5152 9 1 . . 

y = — — 01 f = -((72 cosasm7 -f gi smacos7) , 

ga d-f ' g 

rj 9l92 , 1 . . , . ^ 

Z = 02 f = — [92 sm a cos 7-1-31 cos a sm 7) , 

90- ' 9 



f) 1 

W = — ' ' — — 02 f ^ -(g2sinasin7 — ffi cosacos7) . (B5) 

ga d^ g 



12 

Let us calculate the curvature perturbation for this model. To begin with, we evaluate the perturbation in 7 to the 
second order to obtain 

. ,91,92 ?7i2<^'/>i - mi5(j)2 „ {hlY ^2X - miZ 

g (7 m2Y — miW ' 2 m2Y — miW 

On the basis of these equations, 5N is evaluated to the second order as 



'W5(l)i + Y5(l)2 {9ig2Y {m2Hi - mi5(i)2f 
m2Y — rriiW '^g'^o- {m2Y — miWy 



_ ~yy uyi -r -1 uyJ2 \yiy2) V"'-2"yi — "'-l"y2; ,r^~\ 



Now we can evaluate the quantities of interest. As before, for convenience, we introduce angle (3 as 

51 = 5 cos /? , 32 = .9 sin /3 . (B8) 

Then the curvature perturbation spectrum is 

^5- .^.. r..M. (-) • (B9) 

The spectral index is 



{m2Y~m{WY \2t:J 



The tensor-to-scalar ratio is 



Finally, the non-Gaussianity is 



where we note that 



ns = l-iml+ml). (BIO) 



r-S^-^K^^. (BH) 



Hocai ^ ^ cos^/Jsin /? (miF + maiy)" 
•^^^ 6ct (y2 + vt^2)2(^2r-miVK) ' ^ ^ 



Y^ + W^ = cos^ /3 cos^ 7 + sin^ /3 sin^ 7 . (B13) 



Here, it is worthwhile to mention that the spectral index depends only on rrii and m2. 

To enable a direct comparison with the model discussed in the main text, let us consider the case of VF = for the 
present model as well. In this case, we have 

^' = ^ if)' • 



8m^ 



2 I 



Piocai 55 cos2/3sin 13 



•'^^ 6cr (cos2 /? cos2 7 + sin^ p sin2 7)3/2 ^2 ' ^^^^' 

We see that a large mass ratio, rrii ^ 7712, is necessary in order to realize a large /jvff'. However, because r in the 
present case is determined only by the smaller mass, r = Sm^, it is difficult to realize both large r and large Yn'j^^ 
. This is in contrast to the case we discussed in the main text, for which it was possible to make both values large 
enough to be detectable in the very near future. 
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